Ultrarelativistic black hole formation 
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We study the head-on collision of fluid particles well within the kinetic energy dominated regime 
(7 = 8 to 12) by numerically solving the Einstein- hydrodynamic equations. We find that the 
threshold for black hole formation is lower (by a factor of a few) than simple hoop conjecture 
estimates, and, moreover, near this threshold two distinct apparent horizons first form post-collision 
and then merge. We argue that this can be understood in terms of a gravitational focusing effect. 
The gravitational radiation reaches luminosities of 0.014 c^/C, carrying 16 ±2% of the total energy. 
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Introduction. — An important topic in high energy 
physics that remains poorly understand is the dynamics 
and outcome of super-Planck scale particle collisions. Ac- 
cording to general relativity, kinetic energy, like all forms 
of energy, gravitates. This implies that at sufficiently 
high center of mass energies E the gravitational force 
will eventually dominate any interaction. Suppose one 
can localize the particles' wavefunctions at the moment 
of interaction to be within a sphere of radius R] then, 
according to Thome's hoop conjecture p| (see also [2I- 
3]), if R is less than the corresponding Schwarzschild 
radius Rs = 2GE/c^, the gravitational interaction will 
be so strong that a black hole (BH) will form. For par- 
ticles satisfying the de Broglie relationship the thresh- 
old for BH formation occurs at Planck energies. There 
has been much interest in the past decade over the pos- 
sible relevance of this to proton collisions at the Large 
Hadron Collider [HI, @ and cosmic ray collisions with the 
earth's atmosphere [Tf, spurred by theories of quantum 
gravity with small or warped extra dimensions |8MlO| that 
present the possibility of a true Planck scale within reach 
of these processes. To date no evidence for BH forma- 
tion has been found [ll|, [III , though since the theories do 
not make firm predictions for what the true Planck scale 
is, the high energy scattering problem is worthy (beyond 
intrinsic theoretical interest) of further study. 

Here we explore the purely classical gravitational 
properties of head-on ultrarelativistic collisions (in 4- 
dimensional asymptotically flat spacetime). This ostensi- 
bly gives the leading order description of the process for 
energies sufficiently above the Planck scale, as all non- 
gravitational interactions will be hidden behind the event 
horizon, implying that the particular model for the par- 
ticles is irrelevant. However, part of the motivation for 
this study is to test this notion, and begin to investigate 
how it breaks down approaching the threshold of BH for- 
mation (though again only at the classical level). 

There have been several studies of ultrarelativistic col- 
lisions using BHs as model particles. Penrose [HI first 
considered the head-on collision of two Aichelburg-Sexl 
metrics [l^ , each representing the boost 7 ^ 00 limit of 
the Schwarzschild metric (letting the mass M go to zero 
such that the energy E = 7M is fixed, and note unless 
otherwise indicated we use geometric units G = c = 1). 



Though the spacetime to the causal future is unknown, 
a trapped surface is present at the moment of collision, 
giving an upper bound of 29% for the radiated energy. 
Perturbative methods [H, allowed a direct calcula- 
tion, estimating 16.4% energy emitted. In [l^, head-on 
collisions up to 7 ~ 3 were studied using numerical so- 
lution of the field equations; extrapolating the results to 
7 ^ 00 gave a value of 14 ± 3%. We briefly mention that 
studies of BH colHsions for general impact parameters 
using the trapped surface method for infinite boosts [l8| , 
and numerical simulations of finite boosts lH, 2o| show 
that considerably more energy can be radiated then. 

However, head-on BH collisions at any speed will nec- 
essarily produce a larger BH, and therefore are not suit- 
able for studying the threshold or dynamics of BH for- 
mation, nor whether BH formation is the generic out- 
come regardless of the nature of the colliding particle. 
In [21] a first attempt to address some of these ques- 
tions was made, where the ultrarelativistic collision of 
boson stars (solitonic solutions of a minimally coupled 
complex scalar field) was studied numerically. It was 
found for boson stars with compactness 2M/r ^ 1/20 
that a BH forms for boosts greater than 7 ^ 2.9, roughly 
one- third the value 7h = 10 predicted by applying the 
hoop conjecture at the time of collision. Whether the 
threshold is generically such a factor smaller than the 
hoop conjecture estimate was unclear, first because only 
a single matter model was considered, but also because 
though for 7 = 2.9 there is almost twice as much ki- 
netic as rest mass energy in the spacetime, this may not 
be high enough for the matter dynamics to be irrelevant. 
Furthermore, due to difficulties disentangling gauge from 
gravitational wave (GW) dynamics, no estimates of the 
radiated energy were made. 

In this letter we also study black hole formation 
in head-on particle collisions. However, we use self- 
gravitating, perfect fluid "stars" as the model particles. 
To begin with, this allows us to further test the gener- 
ality of the above arguments in a case where the gravi- 
tational attraction would be opposed by the tendency of 
the fluid to become highly pressurized on collision and 
disperse. Second, the nature of fluid stars, not having 
small-scale internal oscillations as boson stars, as well as 
a new method for constructing initial data j2j|, permits 
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us to explore significantly higher boost collisions where 
the ratio of kinetic to rest mass energy is of order 10:1. 
An independent work with the same matter model used 
here was recently presented in [23|, though as with 21 1 
it focuses on regimes where this ratio is at most ^2:1. 

We find that BHs are formed above a critical boost 7c 
that is a factor of a few less than the hoop conjecture es- 
timate. A new phenomenon we present here is, for boosts 
slightly above 7c, we observe two separate apparent hori- 
zons (AHs) form shortly after the collision, which some 
time later are encompassed by a single horizon that rings 
down to a Schwarzschild BH. We argue that this can be 
qualitatively understood as due to the strong focusing of 
the fiuid elements of one star by the boosted spacetime of 
the other, and vice versa, using a geodesic model similar 
to that in for BH formation in the scattering prob- 
lem. We also study the GWs emitted in this regime for 
the first time and find that for the 7 = 10 BH forming 
case 16 ± 2% of the energy of the spacetime is radiated 
(the energy extrapolation described in [l^ suggests this 
should be 94% of the 7 = 00 limit). For sub-threshold 
cases, the strong focusing leads to high fiuid pressures 
that cause the stars to explode outward. 

In what follows, we outline the equations we are solv- 
ing, the numerical methods for doing so, and the setup 
of the initial data. We then present the results of our 
simulations, compare them to geodesic focusing, and end 
with concluding remarks. 

Methodology. — We numerically solve the Einstein 
field equations, in the generalized harmonic formulation, 
coupled to a perfect fiuid using the code described in [25| . 
For simplicity we use the F = 2 equation of state. We 
use a variation of the damped harmonic gauge [2l|, [26| 
that corresponds to equation (A15) in [26| with p = 1/4. 

We take advantage of the axisymmetry of a head-on 
collision to reduce the numerical grid to two dimensions 
and use seven levels of mesh-refinement where the finest 
level covers the equatorial and polar radii of the star 
by approximately 830 and (due to Lorentz contraction) 
830/7 points, respectively. For the 7 = 10 case, to es- 
timate truncation error we also ran simulations with 1.5 
and 2 times the resolution. Unless otherwise stated, re- 
sults from this case are from the high resolution run. 

Initial data is constructed using free data from two 
identical, boosted solutions of the Tolman-Oppenheimer- 
Volkoff equations with a polytropic condition, and then 
solving the constraint equations in the conformal thin- 



sandwich formulation as described in [22|. With this 
method, the "spurious" gravitational radiation is much 
smaller then the physical signal (see Fig. [3j). We choose 
isolated star solutions with compaction 2M^/R^ = 1/40, 
where and are the gravitational mass and ra- 
dius, respectively, of the star in its rest frame. They are 
boosted towards each other with Lorentz factor 7, at an 
initial separation oi d = 534M*. We consider cases with 
7 = 8, 8.5, 9, 9.5, 10 and 12, though most of our detailed 



results are from 7 = 8 and 10. 

Results. — We find that BH formation does occur 
in the ultrarelativistic collision of fiuid particles with the 
aforementioned compaction for boost factors 7 > 8.5±0.5 
(the uncertainty is from the sampling resolution of our 
survey in 7). This is ~ 2.4 times smaller than the hoop 
conjecture threshold of 7h ~ 20. In Fig.[T]we show snap- 
shots of the rest mass density for a subcritical case with 
7 = 8 and for a supercritical case with 7 = 10. In the 
former, after the collision, the matter focuses down into 
two high density regions which then explode outward. In 
the latter, instead of exploding, two identical AHs appear 
surrounding these regions. (It should be noted that the 
existence of the initially disjoint AHs does not preclude 
the possibility of a signal encompassing event horizon.) 
The AHs then fall towards each other with a third, en- 
compassing AH appearing afterwards. 

In Fig. [2] we show the irreducible mass of the AHs, 
proper distance between the smaller AHs, and the ratio 
of the proper equatorial and polar circumferences Ceq/ Cp 
for 7 = 10. The two smaller AHs are born rather prolate 
with Ceq/Cp ~ 0.6. Together they have mass > 0.4M 
where M ^ 27M* is the total spacetime mass, i.e. they 
contain a significant amount of what was originally ki- 
netic energy. When the third encompassing AH appears 
it initially has less irreducible mass (though greater area) 
than the sum of the smaller AHs. It is also extremely dis- 
torted with Ceq/Cp ^0.2 and an equatorial circumfer- 
ence that is less than the smaller AHs, suggesting more of 
a dumbbell shape. This contrasts with what is found in 
ultrarelativistic black hole collisions where Ceq/Cp ~ 1.5 
initially [l^, consistent with a disk shaped AH. 

In Fig. [3] we show the GW power associated with dif- 
ferent spherical harmonics for 7 = 10, and the early part 
of the GW power for 7 = 8. (Because of the symmetries 
here only the even /, m = harmonics are nonzero.) For 
7 = 10, 16 ± 2% of the initial spacetime energy is ra- 
diated as GWs, with a peak luminosity of 0.0137 ± 1% 
(the error bars include estimates of the truncation error 
and finite radius extraction effects). The mass of the fi- 
nal BH is ^ 0.72M, suggesting the remaining 12% of the 
energy is carried off by the ^ 32% of the initial rest mass 
that remains outside the final BH by the end of the sim- 
ulation. Measuring the contributions to the total energy 
from higher / modes relative to the / = 2 component we 
get that E^/E2 = 0.19 ± 0.01, Eq/E2 = 0.073 ± 0.001, 
and Es/E2 = 0.040 ± 0.002. The substantial amount of 
energy in higher modes is consistent with results from 
ultrarelativistic BH collisions. Also, the zero-frequency 
limit combined with an /-dependent frequency cutoff 
ujc = //(3a/3M) set by BH quasinormal frequencies (27| 
predicts corresponding values of 0.22, 0.09, and 0.05. 

For 7 = 8 we can only extract the GW signal before 
the fiuid outfiow crosses the extraction sphere. Before 
this time, the GW signal looks qualitatively similar to the 
7 = 10 case (see the inset of Fig. [3] for the corresponding 
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FIG. 1. Snapshots of rest mass density on a logarithmic scale from 10~^ to 10^ times the initial maximum density, for 
simulations with 7 = 8 (top) and 7 = 10 (bottom) at times (left to right) t = 0, the initial time; t = 300M*, shortly after 
collision; t = 375M*, after the appearance of the smaller AHs in the 7 = 10 case; t = 424M*, after the appearance the third, 
encompassing AH (white outline) in the 7 = 10 case; and t — 700M*. The black regions are best-fit ellipses to the AHs. 




FIG. 2. Apparent horizon properties for the 7 = 10 case. 
The AHs labeled 1 and 2 are the two identical (when mirrored 
about the collision plane) ones that appear first and are later 
encompassed by AH 3. We show (top to bottom): the irre- 
ducible masses, the proper distance dp between (and exterior 
to) AHs 1 and 2 measured along the collision axis, and the 
ratio of the equatorial to polar circumferences. 

power) and contains 10% of the energy of the spacetime. 

Cases with 7 = 9.5, 9.0, and 8.5 also first form two dis- 
joint AHs with increasing initial separation the smaller 
the boost. However, we were unable to follow these 
cases through merger before numerical instabilities set 
in on the excision surface. The reason, we believe, is the 
smaller boosts form more distorted AH shapes, and our 
current approach of excising the best-fit ellipse to the AH 
shape is inadequate. We have also been unable to obtain 
robust results for significantly higher Lorentz factors due 
to high frequency numerical instabilities that develop at 
the surface of the boosted stars; however it seems that 



FIG. 3. Total and spin- weight 2 spherical harmonic decom- 
position of power in GWs from the 7 = 10 case (with 7 = 8 
in the inset). The small feature at (t — r) < is an artifact of 
the initial data and is not include in the estimate of energy. 

the third AH appears at nearly the same time in this 
gauge as the first two AHs at 7 12 for the stars con- 
sidered here, and for larger boosts we expect a single AH 
to form at collision. 

Geodesic focusing. — To illustrate the manner in 
which a boosted star may act like a gravitational lens 
and, during collision, focus the matter of the other star, 
we consider a simplified scenario in a spacetime consisting 
of a single boosted star. We follow a set of geodesies com- 
ing from the opposite direction with the same Lorentz 
factor, initially distributed to fill out the volume of what 
would have been the other boosted star (i.e., we replace 
the second star with a set of tracer particles). These 
geodesies are shown in Fig. |4] for 7 = 10 with the same 
compaction star (2M^/R^ = 1/40) described above. As 
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FIG. 4. Focusing of a set of geodesies in a boosted star 
spacetime with 7 = 10 and 2M*/i?* = 1/40. Shown are the 
x-coordinate relative to the center of the boosted star which 
is at X* = vt (top) and the ^-coordinate (perpendicular to the 
boost axis) of the geodesies as a function of coordinate time. 



these geodesies pass through the boosted star they be- 
come focused in the direction orthogonal to the boost axis 
while spreading out along the boost axis. The greatest 
focusing (i.e. when the separation between the geodesies 
in either direction is smallest) occurs at a distance of ap- 
proximately 1.5R^ away from the star and reduces the 
radius by a factor of ^ 4 (roughly consistent, with the 
caveat of coordinate differences, with what we observe in 
the full problem — see Figs. [1] and [2]). Contraction by this 
factor is sufficient to get to the BH formation threshold 
if we assume that this focusing also converts sufficient 
translational energy to radial inflow that it is valid to 
apply the hoop conjecture to the total energy of the star 
in this frame. Evidence that this assumption may be rea- 
sonable comes from the temporary slowdown of the trans- 
lational velocity of the geodesies seen in Fig. |4] (though 
somewhat before maximum compression), and from the 
full simulations where in the 7 = 10 case the two AHs 
move toward each other after formation, and in the 7 = 8 
case post-compression the fluid flow is largely radial. In 
the ultrarelativistic limit, this geodesic focusing factor is 
mainly a function of the ratio jM^/R^^ and similar re- 
sults are obtained for larger boosts with correspondingly 
less compact stars. This simplistic treatment of course 
ignores the effects of pressure as well as the non-linear 
gravitational interaction. 

Conclusions. — In this study we considered the head- 
on collision of self-gravitating fluid stars in the regime 
where the ratio of kinetic to rest mass energy in the 
spacetime is ~10:1. We find above a critical boost 
7c = 8.5 ± 0.5 that BHs do form. The dynamics of the 



solution, and a simple geodesic model similar to [2J, sug- 
gest that near threshold the strong focusing nature of the 
spacetime sourced by one boosted star on the other, and 
vice versa, causes the energy to be concentrated post- 
collision around two focal points on axis. In the sub- 
critical case, the material explodes outward from these 
points, consistent with 



2IL |23|; however just supercriti- 



cal we find two distinct AHs that initially form around 
the focal points. This focusing also offers an intuitive ex- 
planation for why the threshold in cases studied to date is 
systematically less than hoop conjecture estimates (here 
7c/ Ih ^ 0.4, with the boson star collisions jc/jh 0-3 
for 7c ~ 2.9 [21], and similar factors were found in |28l430| 
for the scattering problem using a perturbative string- 
inspired model). 

For the 7 = 10 supercritical case, we find 16 ± 2% 
of the total energy is radiated gravitationally, consis- 
tent with results extrapolated from 7 3 BH col- 
lisions [17], and perturbative calculations of the infi- 
nite boost limit [HI, |l6|. Moreover, the leading order 
spherical-harmonic multipole structure of the waves is 
consistent with point-particle approximations and the 
BH case [l^, both super and subcritical, in the latter 
prior to obscuration of the waveform when the matter 
outflow reaches the extraction sphere. 

This suggests three different regimes in the head-on 
collision of ultra-relativistic, non-singular model particles 
in general relativity, for sources that have sufficiently low 
compactness such that 7c ^ 1. For 7 <C 7c gravity plays 
little role, and the dynamics is governed by that of the 
matter; for 7 7c we expect universal behavior, i.e. any 
particle model will give the same quantitative spacetime 
dynamics; however in the intermediate regime 7 ~ 7c 
both gravitational and matter dynamics will be impor- 
tant. Ignoring quantum effects and studying the nature 
of super-Planck scale particles collisions using general rel- 
ativity is arguably robust only when 7 ^ 7c, though 
perhaps some insights can still be drawn from classical 
general relativity in the intermediate regime. 

The intermediate regime includes the threshold of BH 
formation and corresponding, mditteT- dependent critical 
phenomena [3l|. We conjecture approaching 7c may 
generically result in two critical solutions unfolding post- 
collision about the geodesic focal points of the two collid- 
ing particles (we speculate the reason why two distinct 
AHs were not seen in [2l|, |23| is the compactness is not 
sufficiently low to have 7c ^ 1.) For F = 2 fluid stars it 
would be interesting to see whether the critical solution 
is the type I unstable star-like solution found for lower 
7c 's [23], or the type II self-similar solution arising in the 
kinetic energy dominated regime [32[. It would also be 
interesting to explore collisions with nonzero impact pa- 
rameters. This would allow a better comparison to BH 
collisions, which do not have a threshold for BH forma- 
tion, but do have two distinct end states as a function of 
impact parameter: a large BH or two unbound BHs. 
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